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ABSTRACT 

Using recently calculated analytic models for the thermal structure of ultramagnetized 
neutron stars, we estimate the thermal fluxes from young {t ~ 1000 yr) ultramagnetized 
{B ^ 10^*^ G) cooling neutron stars. We find that the pulsed X-ray emission from objects 
such as IE 1841-045 and IE 2259-1-586 as well as many soft-gamma repeaters can be explained 
by photon cooling if the neutron star possesses a thin insulating envelope of matter of low 
atomic weight at densities p < 10 — 10* g/cm'^. The total mass of this insulating layer is 
M - 10-" - lO'^Mp,. 



Subject headings: stars: neutron — stars: magnetic fields — radiative transfer — X-rays: stars 



Introduction 



In recent years, several "breaking" (Mereghetti & Stella 1995) or "anomalous" (Van Paradijs, Taam & 
van den Heuvel 1995 ) x-ray pulsars have been discovered ( Vasisht fc Gotthelf 1997 , Corbet et al. 1995 ). 



These objects typically have pulsed X-ray emission with steadily increasing periods 10 s, X-ray 
luminosities ~ 10'^^ — 10'^^ erg/s, soft spectra, and no detected companions or accretion disks. Furthermore, 
they are typically observed through hydrogen column densities ~ 10^^ cm~^ indicating that they are 



not common. Vasisht fc Gotthelf (1997)| describe several of these sources and present observations of IE 
1841-045, whose properties are characteristic of this class of objects. 



Specifically, [Vasisht fc Gotthelf (1997)| use archival observations of the supernova remnant Kes 73 
obtained with the ASCA and ROSAT satellites. They find that the x-ray source in the center of the SNR, 
IE 1841-045, had a period of 11.766684 s from the ASCA data taken in 1993 October. The ROSAT data of 
1992 March is best fitted with a period of 11.7645 s, yielding a period derivative of P ~ 4.73 x lO^^^s s^^ 
and a characteristic spin-down age of 4,000 yr - close to the estimated age for Kes 73 of 2,000 yr. Using these 
values and assuming that magnetic dipole radiation dominates the spin-down, they estimate the dipolar 
field strength of the neutron star to be ~ 10^^ G, well above the quantum critical field, Bcr ~ 4.4 x 10^'^ G. 
Other anomalous x-ray pulsars (AXPs) generally have small ages and long periods, leading one to derive 
similar field strengths. 

With x-ray luminosities Lx ~ 10^^ — 10'^^ erg/s, AXPs are undcrluminous relative to accretion powered 



X-ray pulsars and are generally isolated. For IE 1841-045, Vasisht fc Gotthelf (1997) estimate a spin-down 
power of 10'^^ erg/s which falls short of the observed luminosity. They also argue that although IE 1841-045 
has a period near the equilibrium spin period for a young pulsar with B ^ 10^^ G and Lx ^ 10"^^ erg/s. 
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only an unlikely evolutionary process could spin down the neutron star to this rate within the 2,000 year 
age of Kes 73. They suggest that IE 1841-045 may be powered by magnetic field decay in a dipolar field of 
strength B ^ 10^^ G (Thompson fc Duncan 1996, Goldreich fc Reisenegger 1992). 



In this Letter, we propose a natural explanation for the observed X-ray emission from AXPs. Neutron 
stars with ages ~ 1,000 yr and magnetic fields B ^ 10^^ G have thermal emission in the X-ray-band with 
total luminosities ~ 10'^^ erg/s, if their surface layers consist of light-weight material, such as hydrogen 
and helium. In previous papers, we have developed an analytic model for ultramagnetized neutron star 



envelopes ( Heyl fc Hernquist 1997a ) and calculated the emission through iron envelopes and showed how 
a strong magnetic field affects neutron-star cooling (Heyl & Hernquist 1997b). Here, we will examine 



the properties of ultramagnetized hydrogen and helium envelopes and draw parallels with the observed 
properties of AXPs. 



2. Model Envelopes 

I 1 

Heyl & Hernquist (1997a) have developed analytic models for ultramagnetized neutron star envelopes 

and find that the flux transmitted through the envelope is simply related to the direction and strength of 
the magnetic field and to the core temperature (Tc). Specifically, these models apply only below densities 
Pmax at which the ground Landau level becomes filled; i.e. 2.2 x 10* — 7.1 x 10^ g/cm'^ for polar fields, 
B = 10^5 - 10^6 G. 

Here we consider hydrogen and helium envelopes with B ^ 10^^ G and core temperatures expected for 
cooling by the modified URCA process after ~ 1,000 yr. Although our models do not extend to the high 
densities traditionally associated with the isothermal core of a neutron star {p ^ 10^° g/cm^^), for B ^ 10^^ 
G and Tc ~ 10* — lO* '"^ K, the envelope is nearly isothermal at /Omax- We denote the temperature at this 
density by Tmax and for this analysis take it to be equal to the core temperature. 

The envelope models are calculated using a plane-parallel, Newtonian approximation. In this approach, 
the core temperature is a function of F/gs, B and ij: (the angle between the radial and field directions). 
F is the transmitted heat fiux, is the surface gravity, and all of these values are taken to be in the 
frame of the neutron star surface. For such strong fields, the models have a simple dependence on the 
angle "0, *-e- F/gs oc cos^ ?A (Heyl fc Hernquist 1997a, Bhibanov et al. 1995, Shibanov fc Yakovlev 1996| ) and 



furthermore the flux for a fixed core temperature is approximately proportional to B^-^. With these two 
facts, we find that the average flux over the surface of a neutron star with a dipole field configuration is 
0.4765 times its peak value at the magnetic poles. 

For a given core temperature, magnetized envelopes transmit more heat than unmagnetized ones 



{e.g. Hernquist 1985). Furthermore, as Hernquist & Applegate (1984) found, the relationship between the 



core temperature and the transmitted flux is strongly sensitive to the composition of the degenerate portion 
of the envelope, specifically iron insulates the core much more effectively that hydrogen or helium. Because 
of the intense magnetic field, the luminosities for a given core temperature are even larger than those found 



by Potekhin, Chabrier fc Yakovlev (1997) 
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3. Luminosity Evolution 

We have calculated the expected luminosity of an ultramagnetized neutron star as a function of time. 



We assume that the core of the neutron star cools by only the modified URCA process {e.g. Shapiro & 

-1/3 f 



Teukolsky 1983) 



At(URCA) ^ 1 H- Tjif) 1 



yr, (1) 



where Tx — K and At is an estimate of the age of the neutron star. For simplicity, we will assume 

that the initial temperature is much larger than the current temperature Tc{f). 

If one ignores gravitational redshift effects which depend on the radius of the neutron star, the photon 
luminosity is given by 

= AttGM- « 9.5 X 1032^^i!Lerg s'^ (2) 

5s 5s, 14 Mq 



We will take M = 1.4M, 



0- 

Figure ^ depicts the photon luminosity and mean effective temperature as a function of time for several 
field strengths with iron, helium, and hydrogen envelopes. We see that the photon luminosity increases with 
increasing magnetic field, but the composition of the envelope is a more important effect. The luminosity 
through a hydrogen envelope is ten times larger than through iron and 1.6 times larger than through 
helium. The relationship between core temperature and flux is most sensitive to the thermal conductivity 



in the liquid portion of the degenerate envelope (Gudmundsson, Pethick & Epstein 1982). In this zone, the 
conductivity is proportional to Z^^; consequently, low Z envelopes conduct heat more readily. 

To lowest order, for a given magnetic field strength and core temperature, the transmitted flux is 
proportional to Z~'^l^ . The dependence on Zj A is weaker. The flux dips sharply for hydrogen and helium 
with B = 10^^ G, when the material near the degenerate-non-degenerate interface begins to solidify. At a 
given density the conductivity due to degenerate electrons is smaller in the solid state than in the liquid 
phase. 

We note that our assumption that the magnetic field has no explicit influence on the phase of the 
matter is problematic for these intense fields. For example, the magnetic field may well alter the state 
of the material in the solid portion of the crust, particularly when the magnetic stress exceeds the yield 
stress of the lattice. The shear modulus of the lattice is ^ (Ze)^n^''^, where nz is the ion number density 



(Ruderman 1972). Assuming that the lattice breaks when the strain angle Xstrain 10 , magnetic stresses 



dominate when {e.g. Blandford & Hernquist 1982; Blandford, Applegate fc Hernquist 1983| ) 



/ \ 2/3 

B > 1.2 X 10^^^1-2/3 ( ) (3) 

which is satisfied throughout the crust for fields B ;> lO^"' G. 

Even for the weakest field considered, B = 10^^ G, the photon luminosity through a hydrogen envelope 
is 1.5 X 10'^^ erg/s for an age of 2,000 yr. This is comparable to both the pulsed X-ray pulsed luminosity 
of IE 1841-045 of Lx ~ 5 x lO'^'^df erg/s {d-j is the distance to the source divided by 7 kpc) and the total 
X-ray luminosity of 3.5 x lO'^^d^ erg/s, given observational uncertainty. 

The photon luminosity is nearly large enough to account for the luminosity in the quiescent state of 
the soft gamma repeater (SGR 0526-66). Several authors have advanced the view that strongly magnetized 
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neutron star power SGRs. For example, rhompson & Duncan (1995) propose that soft gamma repeaters 
(SGRs) are powered by magnetic reconnection events near the surfaces of ultramagnetized neutron stars. 
Ulmer (1994) finds that a strong magnetic field can explain the super-Eddington radiation transfer in 



SGRs. Rothschild, Kulkarni fc Lingcnfclter (1994)| estimate the luminosity of SGR 0526-66 in the quiescent 
state to be approximately 7 x 10'^^ erg/sec. Since SGR 0526-66 is located in a supernova remnant, they 
can also estimate the age of the source to be approximately 5,000 years. For this age and a magnetic field 
of B = 10^^-^G, we find a photon luminosity of 1.5 x lO'^^ erg/sec, not far short of the SGR's quiescent 
luminosity. 



4. Discussion 

To determine the viability of our model for AXP emission, we must ascertain whether or not neutron 
stars with B ~ 10^^ G are likely to have an insulating envelope of the required mass in low atomic weight 
matter. From numerical experimentation, we have found that the relationship between the core temperature 
and the effective temperature is most sensitive to the composition in the degenerate, liquid portion of the 
envelope; i.e. the "sensitivity" region found by [Gudmundsson, Pethick fc Epstein (1982) . The density in 
this zone is approximately in the range 10^ — 10^ g/cm'^. 

We can appeal to the equation of hydrostatic equilibrium to calculate the total mass up to a certain 
density. If we assume that the envelope is thin, we have 

AM . i^P, (4) 

where R is the radius of the neutron star and P is the pressure at the given depth. In the region of interest, 
the pressure is supplied by degenerate electrons occupying the ground Landau level. In a strong magnetic 
field, electrons become relativistic at 

p « 3.3 X 10"Bl5Aieg/cm^ (5) 
where B15 = B/10^^ G and /ig = A/Z is the mean atomic weight per electron {e.g. Heyl fc Hernquist 1997a ). 
For nonrelativistic electrons the equation of state is 

Pe 1.5 X 10^°Pi5VrV6dynes/cm^ for pe < 33Pi5Afe, (6) 

where pe — g/cm"^. 

For densities pe S> 33Pi5/ie, the electrons are relativistic, and we obtain the equation of state, 

Pe = 7.5 X lO^^PfgVrVedynes/cm^ (7) 
Therefore, depending on the maximum density of the insulating layer, we obtain for its total mass, 

AM = 7.1 X (^-^^ B-^p-'^pIMq (8) 

in the non-relativistic limit where M is the mass of the neutron star and Rq = R/10^ cm. In the relativistic 
limit, we obtain 

AM = 3.5 X 10-i3i?4 b^^p-'pIMq. (9) 
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At a minimum, the \ow-Z insulating layer must extend into the degenerate portion of the envelope where 
pe ~ 10- The non-relativistic estimate for the mass of this layer yields ~ 1O~^^M0 for B = 10^^ G. If the 
layer extends to the density at which the first Landau level fills (pe ~ 200) well into the relativistic regime, 
AM - 1O~*^M0 for B = lO^^ G. 

The high temperatures accompanying core collapse are generally thought to process the material that 
will comprise the neutron star envelope to nuclear statistical equilibrium; therefore, the envelope is expected 



to consist of ^^Fe and heavier nuclei at higher densities {e.g. Shapiro fc Tcukolsky 1983 ). A hydrogen or 
helium envelope must contain material accreted by the neutron star from the interstellar medium or from 
the "detritus" of the supernova. Even with an intense magnetic field to dramatically increase the Alfven 
radius, a neutron star would require ^ 10^ yr to accrete ^ XQ^^^Mq from the ISM, unless the neutron star 
were traveling through an unusually dense region. For example. 



M r. 10-^-ny'Rt^'BU\y' ^-^j Mo/yr, (10) 

where is the number density of the ISM in units of lO'^cm^^, — u/10^ cm/s, and v is the velocity of 
the neutron star. In such a dense medium, the neutron star could accrete the minimum required IO^^^Mq 
within w 1,000 years. If the nascent neutron star received a kick during the supernova explosion, it could 



M 
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accrete material from the high density remnant of the progenitor star. Duncan & Thompson (1992) argue 
that a strong magnetic field B ~ 10"'^^ G can cause the neutrino emission from the nascent neutron star 
to be anisotropic; consequently, the neutron star would recoil at v ^ 100 km/s relative to center of the 
explosion. 

I I 

Chevalier (1989) modeled the accretion by a nascent neutron star or black hole from the exploded 
envelope of SN 1987A. He found that ~ O.IMq will fall back within several days of the explosion. For a 
typical type II supernova which has an extended envelope, the total accreted mass may be a factor of 100 
smaller, yielding a total of lO^'^M© of material to comprise the insulating layer. Although much of this 
material will have been processed during the evolution of the star, the X-ray and gamma-ray light curves of 
SN 1987A as well as numerical simulations indicate that heavy and light elements mix during the explosion 



(McCray 1993 and references therein); any unprocessed remnants of the progenitor star's hydrogen envelope 



would quickly float to the surface (Lai & Salpeter 1997) and form a hydrogen envelope on the neutron 



star. Even in an intense magnetic field, pycnonuclear fusion reactions do not proceed quickly enough to 



substantially process the accreted material within 1,000 yr ( Heyl fc Hernquist 1996 ; Lai fc Salpeter 1996 ). 
If insufficient hydrogen fell back, the photon luminosity for a helium envelope is only one third lower than 
that of a hydrogen envelope; consequently, even a Type lb supernova could result in a cooling-powered 
neutron star with sufficient luminosity to be observed as an AXP. 
I I 

Zavlin, Pavlov & Shibanov (1996) calculate several model atmospheres for cooling neutron stars 

including hydrogen and helium for weak magnetic fields {B ~ 10^ — 10^° G), and compare the emergent 
spectra with a blackbody distribution. Because the opacity of a hydrogen or helium atmosphere drops 
quickly with photon energy, high energy photons originate from deeper and hotter layers of the atmosphere; 
consequently, hydrogen and helium atmospheres have spectra which peak at a higher frequency and have 



stronger Wien tails than a blackbody spectrum with the same effective temperature (Romani 1987). Pavlov 



et al. (1996) convolve the spectra from neutron star atmospheres consisting of hydrogen with the ROSAT 
response function and ffnd that the fitted values of T^b are 4.17 times higher than the model value. The 
observed best-fit blackbody temperature for the spectrum of IE 2259-1-586 is 0.45 keV. This is about 
twice the effective temperature of our model envelopes. The processing of the radiation by the hydrogen 
atmosphere may account for this shift. 
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Furthermore, in the energy range of 1-10 keV, hydrogen spectra tend to be less steep {Si, oc 
than the equivalent blackbody spectra. Although they are not as shallow as the non-thermal tail observed 
in IE 2259+586, S^, cx mechanisms in the magnetosphere probably contribute to the high-energy 
emission. 



An intense magnetic field induces a temperature variation across the surface of the neutron star (Hcyl 



fe Hernquist 1997a). This gradient combined with the limb darkening manifest in models of magnetized 



neutron star atmospheres ( Romani 1987 , Zavlin, Pavlov fc Shibanov 1996 ) can naturally explain the 
observed pulsed fraction of Kes 73 of 35% ( Vasisht fc Gotthelf 1997 ). 



5. Conclusions 

We find that young {t ~ 1,000 yr) neutron stars with strong magnetic fields {B ^ 10^^ G) and 
hydrogen or possibly helium envelopes have photon luminosities similar to those observed from anomalous 
X-ray pulsars and soft-gamma repeaters in their quiescent state. The total mass of the insulating layer is 
10^^^ — IO^^A/q. a strongly magnetized neutron star could accrete enough material from the ISM within 
1,000 years only if the ISM is sufficiently dense, n ^ lO'^cm^'^. However, sufficient material is expected to 
fall back onto the neutron star surface following the explosion of massive stars. 
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Fig. 1. — Photon luminosity and mean effective temperature as a function of the age of the neutron star. 
The upper bold curves trace the evolution for hydrogen envelopes, the intermediate curves give the evolution 
for helium envelopes and the lowest light curves follow the cooling through iron envelopes. The values of 
Pmax for hydrogen envelopes are given in the legend. 



